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Abstract 

We discuss the linear response of the Bogolyubov theory of a weakly interacting Bose gas 
below the critical temperature. The specific aspects, due to induced perturbation terms in the 
O ' self-consistent treatment are discussed in detail. We show that Bogolyubov's model having a 

gap-less equilibrium spectrum nevertheless gives rise to an effective susceptibility which does not 
satisfy the super-fluidity criterion of Hohenberg and Martin. 

> 

; 1 Introduction. 

. . The theoretical interest for Bose-Einstein Condensation (BEC) has largely increased in the last years 
i due to spectacular experimental evidence v^ith cold atoms in a magnetic trap j|l| Most of the 
O \ theoretical discussion for weakly interacting bosons is going vv^ithin the frame of the equilibrium Bo- 
^ ■ golyubov model (2]. (See some recent books 111 IH) for the recent status of the theory as well as a 
discussion of this model related to some rigorous results in [B])- This model is actually also a self- 
> . consistent one, being hovv^ever a truncation of the general self-consistent (s.c.) Hartree-Fock (HF) 
^ ! approximation with the purpose to obtain a gap-less quasi-particle spectrum, the former being from 
] the variational point of view thermodynamically better. On the other hand, the evolution in real 
■ time towards equilibrium and the response of condensed bosonic systems to time-dependent exter- 
^ ■ nal perturbations require a consideration of the self-consistent scheme out of equilibrium. Outside 
O . equilibrium the self-consistent averages, including the anomalous ones like the order parameter, are 
time-dependent and give rise to an induced perturbation. This aspect known also in the ordinary 
theory Hartree-Fock theory of Coulomb interacting particles modifies essentialy the theory of linear 
response of a weakly interacting Bose gas. We revisit here the linear response properties of the 
^ ■ Bogolyubov theory . Our analysis shows, that the Martin- Hohenberg [Z| criterion for supra-fluidity 
P . is actually not satisfied in this model. Although the density-density fluctuation at k = diverges 
logarithmically at small frequencies 0, the effective susceptibility remains finite. 



2 The self-consistent Hartree-Fock approximation for weakly 
interacting bosons. 

As it is well-known, the macro-canonical density matrix 
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^ Z 

is the solution of the following variational problem for the free energy: 

G = tnfpTr S^p (^H - fiN + ^ ln(p) 

with p > and Trp = 1 . 

The equilibrium Hartree -Fock approximation is the approximate solution of this problem with 
the variational reservoir 

- r 

where Htest is the most general one-particle test Hamiltonian. 
It is known 0, that if 



Htest - /iA^ = CgAg 



with Aa given operators and the corresponding variational parameters, the minimum of the above 
variational problem is achieved for 

d{H-pN) 

^» = — ^r^\ — where {■■■} =Tr [ptest ■■■) ■ 

If the operators A^ are bilinear and linear combinations of creation and annihilation operators, then 
the (first) Wick theorem for equilibrium averages can be used to express < H — pN > as a function 
of < Aa > for all a 's. This gives rise to the usual HF coefficients. While applying these results to 
bosonic systems one should take into account the existence of anomalous averages. 

Actually the s.c. Hartree-Fock Hamiltonian is determined by the variational principle up to a 
constant ("HF constant") , which may be fixed by the requirement 

{H) = (Htest) ■ 

Here of course the above described Wick 's theorem must be applied also on the left hand side of 
this equation. 

The resulting s.c. Hartree-Fock Hamiltonian may be generalized successfully under the same 
form also to non-equilibrium, however all the non-equilibrium the averages are now time dependent. 
It may be shown, that within this scheme (in the absence of time dependent external potentials), 
although the HF Hamiltonian is time-dependent, the average HF energy is still conserved. 

Let us consider now the Hamiltonian of interacting bosons (in the presence of an external potential 
U{x) ) 



H = J dx |^(x)+ (-^V^ + U{x)^ ij{x) + i y dx'w{x - x')^/'(x)+^(a;')+^(x')^/'(x)| . (1) 

This Hamiltonian conserves the number of bosons, but we are looking for a spontaneous symmetry 
breaking solution. Therefore in applying the HF-scheme we choose for the set Aa the operators 
ip{x)^ip{x'), ip{x)%lj{x')^ ip^x) and their conjugates. 
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It can be shown that, for the case of spontaneous symmetry breaking with anomalous averages 
{i^ix)) and {iIj{x')i(j{x)) , the first Wick theorem needed by the above stated recipe for equilibrium 
gives rise to 



= {tlj{x)+'ilj{x')){ip{x')+'4j{x)) + {ip{x)+tlj{x)){tlj{x')+'ilj{x')) 
+ (^(a;)+^(x')+) {i>ix')^ix)) - 2\{ijj{x)){i>ix'))\^ . 



The last term is the correction to the usual Wick theorem. This result may be obtained also by sep- 
arating 'iIj{x) = ip'{x) + (^(x)) and afterwards applying the usual Wick theorem. After reintroducing 
the old operators one gets the previous result. 

Based on this Wick's theorem one obtains immediately, according to the recipe described above, 
the form of the s.c. HF Hamiltonian Hfjp (the upper index indicates, that it refers to equilibrium) 



Here we took into account also that w{x — x') = w{x' — x) and omitted the HF constant. 

It is important to remark, that the above described procedure differs from a formal shift directly 
in the Hamiltonian followed by applying the operatorial (second) Wick theorem for contractions. 

We shall not discuss here furthermore the existence and subtle definition of the Bose condensation 
in the presence of an external potential (see [H]), but we shall treat later the linear response to a 
time dependent external potential. 

The above described s.c. Hartree-Fock Hamiltonian may be generalized also to non-equilibrium 
situations and in the case of a contact interaction is often called as The Bogolyubov-Popov modeipi 



In the absence of an external potential U{x) we shall rewrite the s.c. HF Hamiltonian in the 
/c-space. One introduces the Fourier transforms 




+ 



+ 



2^p{x)+ip{x){'tp{x'yij{x')) 
{iIj{x)iIj{x'){'i{j{x)'i/j{x'))* + h.c.) 

{A^{x){^j{x)rmx'))\' + h.c.)]} . 



(2) 



cni. 




and one gets due to translational invariance (momentum conservation) 



H 



-0 

HF 




lwo\{bo)m{borbo + h.c.) , 



(3) 
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or 



= Y^i^kbth + kclhb., + c.bt.btn - 2woh{bo)\\{borb, + {bo)b^) (4) 



'V' 

k 



with the notations 



2m V 

p 



Ck = ^^Wk+p{bpb_p) (6) 



V 

p 

3 The s.c. equilibrium HF solution for interacting bosons. 

In the case of a homogenous system in equilibrium one may study in more detail the self-consistency 
equations and therefore the properties of the solution. With the help of the Bogolyubov transforma- 
tion one may bring the Hamiltonian Eq.(m) to the diagonal form 

^HF = 5^ Efcttfc ctfc + const. 

k 

with some new (bosonic) quasiparticle annihilation and creation operators afc,a^. 

1) Keeping only the terms of the HF Hamiltonian containing operators of the non-condensate, 
relevant for the elementary excitations, we have: 

^){btbk + btkb^k) + clbkb.k + Ckbtubt] 

kj^O 

with fi being the chemical potential. 

This Hamiltonian is diagonalized to the form Ylk^o ^kCt^'^k (up to a constant term) by the canon- 
ical transformation ^ to the new boson fields and a_k '■ 

bk = UkOk + VkOitk (7) 

with 

Uk = u_k; Vk = v_k and - = 1 . 

One gets 

Uk = cosh(xfc)e^^^ Vk = - sinh(xfc)e-*<^'=e-^2V-. (g) 

with tfjk defined by 

Ck = \ck\e-'^^' (9) 

and 

tanh(2xfc) = . (10) 

Cfc — 

One has to require of course 
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Cfe 



< 1. 

|efe - /i| 

The phase remains arbitrary and the energies are given by 



Ek = sign{ek - n) \/ {ck - iif - |cfcP . 
For reasons of stability one should have 

efe - /X > 

and therefore 



ek = \/ (efc - - |cfc|2 . (11) 
Both inequalities may be combined to a single one 

Cfe - > \ck\ . (12) 
2) Now one must consider also the part of the Hamiltonian containing only the operators of the 

— * 

condensate {k — 0) 

(eo - fi)b+bo + ^cob+b+ + ^4bobo - 2woVV\P\\P*bo + Pbt) ■ 
Here we introduced the notation 

^^^(U- (13) 
The diagonalization of the bilinear terms is achieved by a similar canonical transformation 

6o = MoOo + v^aQ (14) 



with 



tanh2xo = (15) 

A* - eo 

and the excitation energy is 

Co = V(eo -/i)^ - |coP (16) 

having the similar condition 

Co - > |co| . (17) 
But we have still to shift the operators cxq in order to eliminate the linear terms 

-2wo\/V\P\''{P*uo + PvDoiQ + h.c. 

Then with 

QJo = «o \Puq + Pvq) (18) 

^0 
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we get a complete diagonalization in the form eoQ;,|ao (up to a constant term). 

3) Now we must proceed to discuss the self-consistency requirements stemming from the fact, 
that the system is in macro-canonical equilibrium. Then 

(afc) = and (a+afc) = ^^^^ _ ^ . (19) 



Since {ho) = 'Uo(ao) + vo{ao)* it follows 
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P\mfi-eo)P + c^P*) (20) 



to be a self-consistency equation. 

Since the phase of P may be chosen arbitrarily by a simple phase transformation of the operators 
bo , it is convenient to chose P to be real and positive. Eq. (f2n|l may have a symmetry breaking 
solution (P 7^ 0) satisfying 

^(eo-/i-Co)P^ = l (21) 

^0 

with any real Cq, since for the positivity of eo we had to require already Eq. (fT7ll . 

For eo > 0, -^(oq do) — > for — > oo. Then it follows also in the thermodynamic limit 

l^{b+bo) = \P\' and ^{bobo) = P'. (22) 
The non-condensate parameters in equilibrium are given by 

+ (23) 

and 

f,^ = _|(_i_ + i). ,24) 

Thus the phase — 2?/^^ of Ck differs by vr from the phase of Fk. 

We may go now overall to the thermodynamic limit after separating the p = terms. We have 
in the thermodynamic limit the functions w{k), f{k), F{k), c{k), e{k) and e{k) and therefore using 
the definition of Ck ^(\.^ we get 



or 



c{k) = [ d^pw{k + p)F{p) + w{k)P^ 



with the excitation energies 



:(A;) = ^(e(A;)-^)2-|c(A;)|2 (26) 



Then it follows from the last equation , due to the reality of w{k) also, that c{k) must be real 
(its imaginary part satisfies a homogeneous equation). 
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The parameters of the condensate are determined by 



and the relation Eq. (j2Tjl . 

The chemical potential may be found from the average particle number by 



no + 7^ / d?kf{k) = n . (28) 



(2^) 

(here uq = P"^ ), or 



(27r)3 y V ^(P) (e^^^^P) - 1) 2^ e{p) ^ ^ 



In the case of a contact potential its Fourier transform is constant w{k) = w and as a consequence, 
also c{k) = c with 

c- 

However without a cut-off the integral diverges and therefore c = 

One should mention also, that for T = 0° these HF results coincide with those of the self-consistent 
variational model elaborated by Girardeau and Arnowitt [T2] . 

On the other hand, the Hugenholtz-Pines theorem [T3j (which is a generalization of the Goldstone 
theorem |14| to non- vanishing temperatures ), requires that the spectrum of elementary excitations 
at vanishing momentum starts with vanishing energy. 

Unfortunately, within the s.c. HF scheme the theorem is not satisfied and the elementary exci- 
tations start with a gap, except if |c(0)p = |e(0) — /ip. One may still try to chose a w{q) to satisfy 
this condition P3]. 

The results of this section are well-known in the literature (see [13 lEl for variational approaches 
and other simplified presentations[10j ). Here they are derived using des Cloizeaux's elegant formalism 
for the sake of further reference. 



4 The Bogolyubov model. 

The Bogolyubov theory p] considers a contact interaction between the bosons, but in contrast to 
the full s.c. HF theory, it disregards the fluctuational part of the anomalous propagators. Therefore, 
while obviously from the variational point of view it is poorer, it obeys the Hugenholtz-Pines theorem. 
It may be regarded as the solution of the variational problem with a constraint given by the known 
property of the exact solution. We shall analyze in the next Section, which are the consequences of 
this fact for the linear response from the point of view of the super-fluidity condition of Martin and 
Hohenberg In this Section we just describe the model itself. 



7 



Although the original theory of Bogolyubov [3j was formulated for a spacial homogeneous system, 
we shall define here also its slight generalization in the presence of an external (time dependent) 
potential . 

In the Bogolyubov model one disregards the fiuctuational part of the anomalous propagators 
{{il^{x) — {iIj{x))){'i/j{x) — {ipi^x)))) and considers the s.c. Hamiltonian 

nsit) = J dx l^ip{xy (^-^\/^ + U{x,t) - f?j ipix) 

IX) 

+ — [4?/^(x)+^(x)(V;(a;)+^(x))* 
+ + /i.e.) 

- {Ai,{x){^{x)Ymx))\^ + h.c)]} . (30) 

Then one remains only with three s.c. entities {'tj){x)^tl){x')) , {^{x)) , (■0(x))*. 

In the absence of any potential the equilibrium values of the averages in the homogenous system 
emerge from those of the previous Section after the replacement {ilj{x)ip{x'))o by ('i/')o('0)o = 
Since, we treat here a contact potential we have a unique coefficient c, which in the absence of the 
anomalous propagator is given by 

c = w^{KY = wP^ (31) 

and we choose again P to be real and positive. 
The equilibrium excitation spectrum is 

ek = ^{ek-^iY-\c?. (32) 

where we took into account that e(g) — /i > . 

The analysis of the A; = part of the Hamiltonian gives rise to the s.c. condition for the existence 
of the anomalous solution Eq. (f2T|l . which may be written as 

eo-/i-c=^. (33) 

From the last three equations results that eo = and therefore it follows 

eo - /i = c = wP"^ . (34) 

To conclude, in the Bogolyubov model (in the absence of an external potential) in the thermo- 
dynamic limit 



t{q) = J2wP^^q^ + ] (35) 



2m \ 2m 



Hence the equilibrium excitation spectrum is gap-less and starts linearly with the momentum. 
The average particle number is given by 



^ {2nrJ ""Py e{p) (6/3^^-1)^2^ e(p) J 



n . 



The introduction of a cut-off here is not necessary, since the integral converges for the Bogolyubov 
spectrum, which is linear in p at small momenta and quadratic at large momenta. 
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5 Linear response. 



Although the hnear response to an external potential within the frame the Bogolyubov model was 
already discussed in the literature in some detailjHl we want to reopen the discussion here, since 
important features due to the self-consistency have been ignored until now. The modifications shed 
a new light on the problem of supra-fluidity of the Bogolyubov model. According to Hohenberg and 
Martin [2] a sufficient and necessary condition for super-fiuidity is that the density response to an 
external potential in Fourier-space diverges at vanishing frequency and wave vector. The Hugenholtz- 
Pines theorem [T^ ensures this feature in the frame of the exact hamilton theory. In other words, the 
existence of gap-less quasi-particles in equilibrium is the key for the understanding of super-fiuidity. 
Of course the nature of the singularity may be complicated by the fact, that multi-quasiparticle 
states have also a vanishing threshold. 

The problem is however much more complicated in the case of self-consistent theories, like the HF 
or the Bogolyubov model, since s.c. parameters (averages) are present in the s.c. Hamiltonian and 
they are themselves modified by the external perturbation. In this sense one speaks about an induced 
perturbation. Our next purpose is to analyze this response in the case of Bogolyubov's model. We 
ignore here any dissipation, therefore any statement about singularities of the response theory has 
to be regarded with some caution. The inverse statement is however true: if without dissipation one 
gets no supra-fiuid properties, no such properties will arise due to dissipation either. 

5.1 General formalism 

We consider now the time-dependent self-consistent linear response to a time dependent potential 
U {x, t) starting from the s.c. HF equilibrium in the absence of the external potential. 

We have to take into account however the peculiarities of the "equilibrium" distribution in the 
presence of spontaneous symmetry breaking. The broken symmetry is the particle number conserva- 
tion and therefore, under condensation conditions, H^^ does not commute with the particle number 
operator A^. Therefore the macrocanonical density matrix ^e~^^^HF-t^^)is not a stationary state of 
the system in the usual sense [Hllll]. More precisely, stationarity is reached by the "rotated" density 
matrix: 

which satisfies the Liouville equation 

d ^ 

generated by the "Hamiltonian" 

Indeed, in the spirit of the Van Hove limit of weak coupling to a thermal bath, p{t) approaches, 
for large times, the value ^e~'^^^HF~^'^\ In the absence of the spontaneous symmetry breaking, the 
rotated density matrix is, of course, identical to the density matrix itself. 

We denote 

{{A)),^Trip{t)A) 
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and 

{A)t = Tr{p{t)A) = e-^n--^\{A))u (m = 0, 1, 2, ...) . 

We assumed here that the operator A has the structure . Then m = n — n . Obviously 

the normal entities {n — n) are stationary in the peculiar "equilibrium" state, however the anomalous 
ones {n ^ n) oscillate with a multiple of the chemical potential. 

In the absence of the potential U{x,t) the system is supposed to be in "equilibrium" 



, 1 
p\t=-oo = -^e 



with translational invariant averages 



where {Fk = y '^ky^o e^^'^^^^'^Fk and have the same phase up to a multiple of | and actually both 
may be taken to be real. Further it is assumed, that the self-consistency relations in equilibrium are 
satisfied. 

Again, in order to simplify the calculations, we treat explicitly only the case of a contact particle 
interaction. In this case 



UHpit) = I dx^ip{x)^ {--^V'^ + U{x,t) - h)'iIj{x) 



in 

+ - [^^l;{xy^l^{x){'^{xyi;{x)) 

+ {i;{x)i;{x){ij{x)il}{x)Y + h.c.) 
- {Ai;{x){i^{x)Ymx))\^ + h.c)]] 



According to linear response theory, for any observable A one has 



where 



or in Fourier transforms 



(t) = ei'^"p*Ae~i'^"p^ , 



HF 



^^^'^^^'^ = dte-'{{[A^.^^{t),5nHF{u^)])f. 

Here 5TlHF{t) contains all the terms of first order (direct and induced) of the difference Ti.HFi't) — 
F^HF ■ Thus, taking also into account that (with real P ) 

5 {m{xm\{ij{x)))) = 2PHmx))) + \p^\5mx))Y = p\25mx))) + smxw) 

we have 
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SHHFit) = jdx{ {U{x, t) + 2wS{{iP{x)+iP{x)))t) i;{x)+ilj{x) 

111 

+ ismxr)rMxr + h.c.) 

- 2wP' [{smx)))t + 25mxm) i^{x) + h.c]] . 

In what follows we restrict further the discussion to the Bogolyubov version of the theory, where 
no anomalous correlations are taken into account. Then we have 



Sn^pit) = j dx{{U{x,t)+2w6{{ij{x)+ij{x)))t)ij{x)+ij{x) 

+ '^{2P5mxWMxr + h.c.) 
- 2wP'[{5mx)))t + 25mxm)^l;ix) + h.c.]} . 
In order to simplify the formulae let us define now 

'ilj{x)~^il!{x) \ / 
Pipix) , B = I - 3Ptp{x) 

Pip{x)+ J \ - 3PipIx) + 

and 

5Hf^(t)=^2u) J dxBaix)^5{{Aa{x))) + J dxU{x,t)Bi{x) . 

a 

Then, after Fourier transforming {f{k,uj) — J dx J dte~^^'^^~'^'> f {x , t) ) 

a' 

with 

Xaa'[K,Uj) = — 

Remark: By the interpretation of the susceptibility one has to perform eventually also a shift 
with j-m, (m = 0, 1, 2, ...)in the case of particle number non-conserving terms! 

The relation of the averages of the different entities to the external potential U{k,u!), with 
matrix notations, is given by 

5.2 The structure of x{k,oj) at T = 0°. 

With the notation = (0 < ^ < 1 ) one gets fj, — {2 — ^)nw. We chose 2nw as natural unit for 
energy and k = ^ '^"^^ for momenta. The new variables are : 



/ dt J dxe-<'^--'\{[A^n%,{x,t),B:,{0,0)])r- 



u„ 



1 / 9' + 



\ 2 \y/WTV' 



1 / 9' + le 



2nw 



Since at T ^ 0° we have 



.ep- fJ, 1 







^ (e/^^. -1) 

the self-consistency equation looks at T = 0° in the new units as as: 



(27r)3 n ^ . 

giving rise to the simple dependence between parameter ^ and the density n 

1 6(1-0 



(27r)2 n 

We have to apply the commutations rules: 



[ag, ttp] = 0; [ag, a J] = Sg^p 

[aqapap] = ag{6q^p + (^^y) 

[a^ttg/, ftpttp/] = 0; [agag>,a^ap] = 6p^q>aqap + Sgi^p^agO^ + 5q^pa^,aq' + 5, 
and since at T = 0° 

(0|«g|0) = (0|q;+q;,,|0) = 

the only non-vanishing terms are 

K>«p]|0) = 5g,p; (0|[Q;gQ!g/,Q!+Q!j]|0) = 5p,gV,<Z' + ^P,<?'V,<2 



(36) 



qyCtp Ug' 



and their conjugates. 
Further 



A{k,t) 



Eg«-ik«?«ik-?«9e-<^'=-«-^«)* + ^;q+jk^x*a+a+fc_^e<^'=+«^'')*) + v^P(^^jkajke-'^'=*+^;jfcQ;+^e*^'=* + h.c.{-k)y 
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^(0)+^(0) 
S(0)+ = I i^(0)+2 - 3P^(0)- 
1^(0)' - 3Pij{x) 



2y Z^p,p' ^ Z^p "p Z^p "p 

\ 2 y V^p + 7^ Z]p ^p ~ y 



or 



5(0)- 



F Ep,p'(^pVa-p"p' + «pV"p "-P') + ^ Ep(^p"p + Vpatp + h.c.) ^ 

w Ep,p' "J "p + ^^-p^^-p'ttp'ttp + ?^;^'*p'"p "p') - Ep(«;"p + ^^;"-p) 

\^ W Ep,p'(«P^'P'"p'"p + t'-pt'-p'ajQ;+ + UpV-p'apap) - 2^ Ep(wp"p + Vpotp) J 



At zero temperature T = 0° we have to average over the ground state (vacuum of quasi-particles) 

Xaa'{k,t) ~ (0|K„^.^(a;,t),£+(0,0)]|0)^^ 

XniK t) ~ P'e-^'''\u, + v*_,\' + ^Yl e-'^''-^'''^^%v;_,{ul_^v, + Vk-^u;) - c.c{-k) 



q-k\'^q\ "'k-q 



Xis{k,t) - J] {e<''+''+'''^'\u,W+kU-q-k - e-<'^-''+'''^'\v,-kWv-q) 
q 

-2P2 ((«, + v*_,)v_ke-'''' - {u*_, + Vk)u^ke'''') 



X2iik, t) ~ {ukiv^k + ul)e-''^' - Vkiu^k + vl)e''>'') 
X22{k,t)^2Pmvk\'e^^'''-\uk\'e-^'''') 
X23{k,t) ~ {vkU^ke'''' - UkV^ke-'''') 
XsiiK t) ~ P' {vKvk + u*_,)e-'^^' - uKuk + v*_,)e'^^') 
Xs2{k,t) ~ 2P' {ulvUe^''^' - vluUe-'^'^') 
Xs3{Kt)^2Pmuk\'e^^''-\vk\'e-^''') 



and 
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n(27r) 



dq 



1 



1 



1 



UJ + + €q^k ^ — — ^q-k 



UJ + Ek UJ - €k^ 
UgVq_k{Uq-kVq + Vq_kUq) 



K 



n{2Ti) 



Xuik, u) = 2^{uk + Vk) [ Uk— Vk—^ — 

uj - ek oJ + ek 



^ — — ^q-k 



Vq-kUqUq_k 



1 



+ eq + tq-k 



UqVg_kVq 



K 



n{2Ti) 



Xnik, uj) = 2^{uk + Vk) ( Vk— — 

UJ - tk 



1 



Uk — 



+ ek 



1 



1 



UqVq.kVq 



^ — ^q — ^q-k UJ + tq + Eg^k 



-Vq^kUqUq^k 



X2iik,uj] 
X22{k,uj) = -2^ (^\vk\'^ 
X23{k,uj) = -2^UkVk 



-^{Uk + Vk) Uk- Vk 



UJ - ek 



uj + ek 
1 



UJ + ek 
1 



\Uk\ 



UJ - ek 
1 



UJ + ek UJ - ek 



Xuik.uj) = -^{uk + Vk) [vk— Uk 



X32{k,uj) = -2^UkVk 



UJ - ek + ek 

1 1 



UJ + ek UJ - ek 



X33ik,uj) = -2^ \uk 



1 



1 



\Vk\ 



UJ + ek UJ - ek. 

As it may be seen, the susceptibilities have poles in the complex u plane at and cuts on the 
real axis from 2ek/2 to oo , respectively from —2ek/2 to -oo . Since ek > 2ek/2, the poles are embedded 
in the cuts. 

The density-density correlation xii{k,uj) of the Bogolyubov model was already discussed at ar- 
bitrary temperatures in Ref-jH]. Our accent is on taking into account the induced perturbation and 
therefore needs the knowledge of the whole matrix x{k,uj). 
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5.3 The density-density response at k = and co 

At k = 0, with the notations 



0. 



CliUJ = 



ra(27r)3 

^3 



dq 



2e„ 



2 2 



n(2ff) 



n(2^) 



UJ 



one gets 



ci(cj) 



25 



i 



M + ii 

,7, "T ,7>2 



(37) 



2€ , e 

,-. "I" ,-,2 



Of course the coefficients of the integrals are not independent entities and are related to the 
parameter ^ by the self-consistency relation Eq. (j36|l . 
Remark, that 

Xii(0, lS) ~ logo; for w ^ 

since the residua of its poles vanish at A; = 0. Thus the naive response theory[S], which ignores 
subtleties of the induced perturbation leads indeed to a singular response, although the singularity 
is only a logarithmic one. The cuts responsible for this behavior are due to excitations of two quasi- 
particles. These are however not present in the near-equilibrium oscillations of the solution of the 
Gross-Pitaevski equation. 

However, taking into account the induced perturbation from Eq. (jrr7|l it follows 



Sii(0,a; 



x(o, 



i-x(o,^) 

Since C\{uj), C2(c<;), 



_ -2^^(2 + 3ci(a;) + €2(0;) + c^{uj)) + ^0^(03(0;) - €2(0;)) + u?Ci{uj) 
~ 2e(-l + 3ci(a;) + 02(0;) + CsH) - ^0^(03(0;) - c^iu)) + - Ci(a;)) ' 

(38) 

€3(0;) behave like logu; for a; — 0, it follows 



lim Sii(0, u) 



-1 . 



(39) 



Thus we may conclude, that although the Bogolyubov model has quasi-particles with vanishing 
energy at g — > 0, the response to an external potential is not at all singular! 

The special value —1 of the effective susceptibility implies that the Fourier transform of the 
effective induced potential for the non-condensate defined as U2ff{x,t) = U{x,t) +2w6{'i/j{x)~^ilj{x))t 
vanishes in the limit — > and a; ^ 0. The factor 2 in this definition stems from the fact, that in 
the case of a contact interaction the direct and exchange terms are identical. On the other hand, the 
effective induced potential for the condensate is different Uljj{x,t) = U2fj-{x,t) — w6\{'ilj{x))t\^ and 
does not vanish in this limit. 

Superfluidity would require diverging density response to a finite external perturbation and this 
does not happen. 
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6 Conclusions 



We have shown, that although the Bogolyubov model leads to a gap-less quasiparticle model, if one 
takes into account the induced perturbation, it leads to no low frequency, low wave- vector singularity 
of the response to a perturbation by an external potential. In this sense it does not fulfill the hopes 
to satisfy the criterion of Martin and Hohenberg for superfluidity. Earlier treatments of the linear 
response in Bose condensed systems either ignored the induced perturbation [H] or did not discuss in 
all detail his consequences, their attention being concentrated on improving the Hartree-Fock theory 
[T7| IT8] . As it is known, the Bogolyubov model is also not the best s.c. approximation from the 
thermodynamical point of view, but still obeys the Huggenholz-Pines theorem, which the thermo- 
dynamically better Hartree-Fock theory does not obey. Seemingly none of the simple self-consistent 
approximations is quite satisfactory and better models have to be developed. 

One of the authors (A.M.B) thanks the Deutsche Forschungsgemeinschaft for the generous support 
allowing his stay at the Frankfurt University. 
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